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One Majorana doublet can be realized at each end of the time-reversal-invariant Majorana 
nanowires. We investigate the Josephson effect in the Majorana-doublet-presented junction mod¬ 
ified by different inter-doublet coupling manners. It is found that when the Majorana doublets 
couple indirectly via a non-magnetic quantum dot, only the normal Josephson effects occur, and 
the fermion parity in the system just affects the current direction and amplitude. However, in the 
odd-parity case, applying finite magnetic field on the quantum dot can induce the appearance of 
the fractional Josephson effect. Next, when the direct and indirect couplings between the Majorana 
doublets coexist, no fractional Josephson effect takes place, regardless of finite magnetic field on the 
quantum dot. Instead, the 7r-period current has an opportunity to appear in some special cases. 
All the results are clarified by analyzing the influence of the fermion occupation in the quantum 
dot on the parity conservation in the whole system. We ascertain that this work will be helpful for 
describing the dot-assisted Josephson effect between the Majorana doublets. 

PACS numbers: 03.75.Lm, 74.45.+C, 74.78.Na, 73.21.-b 


I. INTRODUCTION 

Topological superconductor (TS) has received 
considerable experimental and theoretical atten¬ 
tions because Majorana zero-energy modes appear 
at the ends of the one-dimensional TS which can 
potentially be used for decoherence-free quantum 
computation^""— In comparison with the conven¬ 
tional superconductor, the TS system shows new 
and interesting properties'^ For instance, in the 
proximity-coupled semiconductor-TS devices, the 
Majorana zero modes induce the zero-bias anomaly 
A more compelling TS signature is the unusual 
Josephson current-phase relation. Namely, when the 
normal s-wave superconductor nano-wire is replaced 
by a TS wire with the Majorana zero modes, the 
current-phase relation will be modified to be Ij ^ 
sin ^ and the period of the Josephson current vs (j) 
will be Jtt {(j) is the superconducting phase differ¬ 
ence). This is the so-called the fractional Josephson 
effectSuch a result can be understood in terms 
of fermion parity (FP). If the FP is preserved, there 
will be a protected crossing of the Majorana bound 
states at 0 = TT with perfect population inversion. 
As a result, the system cannot remain in the ground 
state as (j) evolves from 0 to 27r adiabatically^iS 

Recently, a new class of topological superconduc¬ 
tors with time-reversal symmetry, referred to as a 
Dill symmetry-class superconductor and classified 
by the Z 2 topological invariant^""— has attracted 
rapidly growing effortDifferently from chiral 
superconductors, in Dlll-class superconductors, the 
zero modes come in pairs due to Kramers’s theorem. 
Up to now, many schemes have been proposed to re¬ 
alize Z 2 time-reversal-invariant Majorana nanowires 
using the proximity effects of d-wave, p-wave, si- 
wave, or conventional s-wave superconductors 
It was shown that at each end of such a nanowire 
are localized two Majorana bound states that form a 


Kramers doublet and are protected by time-reversal 
symmetry. It is certain that the Majorana doublet 
can induce some new and interesting transport be¬ 
haviors. Some groups have reported that in the 
Josephson junction formed by the Majorana dou¬ 
blet, the Josephson currents show different periods 
in the cases of different FPs^^i However, for com¬ 
pletely describing the transport properties of the 
Majorana doublet, any new proposals are desirable. 

In this work, we aim to investigate the influence 
of an embedded quantum dot (QD) on the current 
properties in the Josephson junction contributed by 
the Majorana doublets. Our motivation is based 
on the following two aspects. First, QD is able to 
accommodate an electron and the average electron 
occupation in one QD can be changed via shifting 
the QD level. Thus, when one QD is introduced 
in the TS junction, the FP can be re-regulated and 
the fractional Josephson current can be modified. 
Moreover, some special QD geometries can induce 
the typical quantum interference mechanisms, e.g., 
the Fano interference^ which are certain to play an 
important role in adjusting the fractional Joseph¬ 
son effect. Second, one QD can mimic a quantum 
impurity in the practical system, which is able to 
provide some useful information for relevant experi¬ 
ments. Our calculations show that when the Majo¬ 
rana doublets couple indirectly via a non-magnetic 
QD, only the normal Josephson effects take place, 
irrelevant to the change of the FP. When a finite 
magnetic field is applied on the QD, the fractional 
Josephson effect comes into being in the odd-FP 
case. On the other hand, when the direct and in¬ 
direct couplings between the Majorana doublets co¬ 
exist, no fractional Josepshon effect occurs despite 
the presence of magnetic field on the QD, and the 
periods of the Josephson currents have opportunities 
to present the 7r-to-27r transition when the QD level 
is shifted. The results in this work will be helpful for 
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FIG. 1: The Josephson junction formed by the direct 
coupling between the Majorana doublets and their indi¬ 
rect coupling via a QD. 

describing the QD-assisted Josephson effect between 
the Majorana doublets. 

II. MODEL 

The Josephson junction that we consider is formed 
by the direct couplings between the Majorana 
nanowires and their indirect couplings via a QD, as 
illustrated in Fig.l. The particle tunneling process 
in this junction can be described by Hamiltonian Ht 
with 

Ht = ^ HaMH to + Hti- (1) 

a=L,R 

HaM denotes the particle motion in the two Ma¬ 
jorana nanowires. With the proximity-induced p- 
wave and s-wave superconducting pairs, the effective 
tight-binding Hamiltonian in the a-th nanowire can 
be written as^ 

HaM 

= '^^tajcl^j^CaJTlcr + + H.C. 

jo- j 

^ ^ap^a,ji^a,j+U ^ 

j 

+ (2) 

j jo 

c^a.jo- ^cx,j(T {o- =t,i or ±1) are the electron cre¬ 
ation and annihilation operators for the j-th site in 
the a-th nanowire, tja is the inter-site hopping en¬ 
ergy and ta,so represents the strength of spin-orbit 
coupling. Aap and denote the energies of the 
p-wave and s-wave superconducting pairings, respec¬ 
tively. Pa is the chemical potential in the a-th 
nanowire. Note that the hopping coefficients and 
the chemical potential are generically reonormalized 
by the proximity effect. The second term Hto de¬ 
notes the direct coupling between the two Majorana 
nanowires, which can be expressed as 

Hto = E + H.C., (3) 

a 

where T is the direct coupling coefficient. Next, Hti 
is to express the indirect coupling between the two 


Majorana nanowires due to the presence of an em¬ 
bedded QD (or a quantum impurity). Its expression 
can be given by 

Hti = T! + R{djd^ + d^d-f) + U nd^Udi 

a 

+ ^ VLdlcL,Na + ^ yndlcRpo- + H.<^A) 

a a 

Here dj. and are the electron creation and annihi¬ 
lation operators in the QD, and eo is the QD level. 
R denotes the strength of an effective magnetic field 
applied on the QD, and U denotes the intradot elec¬ 
tron interaction with Uda = dj.dcr. In addition, Va 
is coupling coefficient between the QD and the a-th 
Majorana nanowire. 

In order to discuss the Josephson effect in this 
junction, we have to deduce an effective Hamilto¬ 
nian that reflects the direct and indirect couplings 
between the Majorana doublets. For this purpose, 
we define the Majorana operators 

7ai = 

j 

7al = 

j 

7a2 = Et^S' 
j 

7a2 = X!+ Mai (5) 

j 

where Caj = SiCajf + ^ 2 Ca,ji is the renormalized 
electron operator at the j-th site in the a-th site 
with Caj = TcajT~^- Using the above formulas, 
we can solve the electron operators in terms of Ma¬ 
jorana and nonzero-energy quasiparticle operators. 
Reexpressing the quasiparticles in terms of electron 
operators, we can interpret crn^ cri, and cri 
by 

cln = Mln*7l2 -'^aLjCLj - T] ^bdi’ 
j j 

Cln = Mln*7l2 - dLjCLj - ^Licfi’ 
j j 

Cri = Mm*7m - '^o-RjCRj ~ E^b^b’ 
j j 

Cri = Mm*7m - '^d-RjCRj - E^b^b’ 
j j 

in which the normalization factor has been ne¬ 
glected. Besides, aaj^ daj and haj, haj are expansion 
coefficients, originated from the quasiparticle opera¬ 
tors other than the corresponding Majorana mode. 
Substituting Eq. into the expression of Ht^ we 
can obtain the low-energy effective Hamiltonian of 
Ht in the case of infinitely-long nanowires, which is 
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divided into two parts. The first part is 

= iFocos ^( 7 L 27 m - 1l21ri) + Ssdld, 

+Unsns - iWLe-^^/^dljL2 + iWLe-^'^^^dljL2 
+Ssdtds + Wndl^fRi + Wndl^Ri + H.c.. (7) 


the relations in Eq.(j6]), we can get the relationship 
that a^jhaj = \ciaj^aj\ and b^jhaj = • 

Accordingly, can be written as 

"Hy = zrIL sin (hL2lL2 - «r IR sin (jr-fRilm 

+(r2Le-*^ + r2fi)44 + H.C., (10) 


The relevant parameters here are defined as fol- 
lows: To = Wr = mdArl> and 

Wr = in which 

Alw = and 

Ai?i = (It is reasonable to suppose 

IMlnI = IAlwI and Besides, in the 

above formula dl = + s^dj), d~ = (—s^dj + 

s^dj), and Sg/s = Sq ^ R with Us = d|ds and 
ris = dtds. 

For the second part, when the highest-order terms 
are neglected, it can be approximated as 

lj!'llilL 2 CZ^aRjCRj + 

3 3 

3 3 

-'^f^LNlL2(Y^dRjCRj + y] b*Rjc\^j) 

3 3 

-Tf'Am 7 m(F ^^jCLj + 4 i 4 i) 

3 3 

~ dajCaj + 

3 3 

- E ^- 4 (E + E ^-. 4 )+( 8 ) 

3 3 

It should be noted that since the s-wave pairing is 
present in the quantum wires, the electrons and 
Ca will form a Cooper pair and condense. This pro¬ 
cess leads to an effective coupling between Majorana 
zero modes localized at the same end and the finite 
coupling between the Kramers doublet in the QD. 
Therefore, up to the second-order perturbation in 
the tunneling process, we can express as 

= ^{'^^fJ-LNdLN'yL 2 lL 2 + V| 44 )^fl 

+ ^(^Vmbm7i?i7m + Vldldl)gL + H.c.,(9) 

where Qa = J2j (^ajdaj J dT{TrCaj{T)Caj{0)) + 

Ey44i/c^^(TVc4(44i(0)) with JdT{---) 
being a time-ordered integral. In the case 
of uniform superconducting pairings in the 
Majorana nanowires, Qa can be further de¬ 
duced as Ga = CLajdaj ^^^(AaT,Aap) ~ 
Ej Kjkj Efc 

eigen-energy of the isolated superconductor With 


in which Til = Bia = 

and r2a = Up to now, we have obtained 

the low-energy effective Hamiltonian of such a struc¬ 
ture. 

The phase difference between the two Majorana 
wires will drive finite Josephson current through 
them, which can be directly evaluated by the fol¬ 
lowing formula 


2e {ORt) 
h dcf) 


(11) 


with (• • •) being the thermal average. It is certain 
that solving the Josephson current is dependent on 
the diagonalization of Rt- 

In the following, we try to diagonalize the Hamil¬ 
tonian. First, by defining 71 = -^{jl2 +7i?i) and 

72 = -^(7i?i+7i2) with 7j = T'yjT~^, we reexpress 
Ht as 


%T = i(ro cos^ +risin(/)) 7 i 72 + 

-i(ro cos ^ - Ti sin 0)7172 + e~sd\ds 

+Un,ns + (r2Le-*^ + r 2 i ?)44 

+ T(_iWie -*^/24 + VFfl 4 ) 7 i 

+ + WRd\)-~i 2 

V2 

+ ^{iWLe-^^/^d\ - WRd\)^i 

+ ^{iWLe-^^/^dl + B^fl4)72 + H.C., (12) 

V2 

where Tiq, is supposed to be Ti. Next, Rt can be 
expressed in the normal fermion representation by 
supposing 7i = (/ + p), 72 = i(/l - /) and ji = 
- /)) 72 = if + P) where p, p and /, / are 
the fermionic creation and annihilation operators. 
Accordingly, the matrix form of Rt can be deduced 
on the basis of InsrisUfrij) where n/ = f and = 

P f. Note that in the system with Majorana bound 
states, only the FP is the good quantum number. 
Thus, we should build the Fock state according to 
the FP. First, in the case of the even FP, the Fock 
state can be written as |Te) =ai| 0000 )-l-a 2 | 0011 )-l- 
a3|0101)+a4|1001)+a5|0110)+a6|1010)+a7|1100) + 
aglllll). As a result, the matrix form of R^^ can 
be written as 
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-A 

0 

-A 

—iA 

-A 

iA 

V 

V 

0 

A 

B 

iB 

-B 

iB 

V 

V 

-A* 

B* 

£g ft 

0 

0 

0 

iB 

—iA 

iA* 

-iB* 

0 

Sg ft 

0 

0 

B 

-A 

-A* 

-B* 

0 

0 

Eg + fi 

0 

-iB 

—iA 

-iA* 

-iB* 

0 

0 

0 

£g fj 

B 

A 

v* 

V* 

-iS* 

s* 

iB* 

s* 

£g Sg U — A 

0 

V* 

V* 

iA^ 

-A^ 

iA* 

A^ 

0 

Eg Ss U “h A 


where A = - Wr)/V2, B = + 62 IOOIO)+ 63 IOIOO)+ 64 IIOOO)+ 65 IOIII)+ 66|1011) + 

hFi?)/A/2, V = r 2 Le*'^ + r 2 i?, A = 2rocos|, and 67 11101) + 63 11110) and the matrix of takes the 

ft = 2risin(/). Next, for the case of the even FP, form as 

the Fock state can be written as |4^o) = ^i|0001) + 


J 


-fi 

0 

B 

-iB 

-A 

iA 

V 

V 

0 

fi 

B 

iB 

A 

iA 

V 

V 

B* 

B* 

Eg A 

0 

0 

0 

iA 

—iA 

iB* 

-iB* 

0 

Sa- A 

0 

0 

-A 

-A 

-A* 

A* 

0 

0 

Eg + A 

0 

-iB 

-iB 

-iA* 

-iA* 

0 

0 

0 

£g A 

B 

-B 

V* 

V* 

-iA^ 

-A* 

iB* 

S* 

Eg Es E U — ft 

0 

v* 

V* 


-A* 

iB* 

-S* 

0 

Eg E Eg E U + U 


For the extreme case of strong magnetic field limit, 
if Eg is in the finite-energy region. Eg will be empty, 
and then only one level contributes to the Josephson 
effects, respectively. Accordingly, in such a case, the 
matrixes of and will be halved, i.e.. 


- -A 

0 

-A 

-A 1 

0 

A 

B 

-B 

-A* 

B* 

Eg ft 

0 

_ -A* 

-B* 

0 

Eg -f ft _ 


(15) 


and 


E 


(o) 

T 


- -ft 0 B -A - 

0 n B A 

S* S* 55-A 0 

_ -A^ A^ 0 55 + A _ 


(16) 


Based on the above analysis, we can understand 
that at the zero-temperature limit, the Josephson 
current in this structure is dependent on the FP, 
i.e.. 


(e/„) ^ ^ 

h d(j) h dcj) 


(17) 


Eqq^^ are the ground-state (GS) energies in the even 
and odd FP cases, respectively. 


III. NUMERICAL RESULTS AND 
DISCUSSIONS 

Following the derivation in the above section, we 
next perform the numerical calculation to discuss 
the detailed properties of the Josephson current 
through such a system. As a typical case, only the 
zero-temperature limit is considered. Before calcula¬ 
tion, we would like to review the Josephson effect in 
the case of Vk =0. In such a case, Et = (Fq cos ^ + 
Fi sin0)(2n/ — 1) + (Fq cos ^ — Fi shi(j)){2nj — 1), 
and \nfnj)f are the eigenstates of Et- The two 
even-FP eigenstates are |00)/ and |11)/, and their 

corresponding GS energies are Eq\ = ^2rocos|. 
Gontrarily, the odd-FP eigenstates are |10)/ and 

|01)/ with the GS energies Eq\ = ±2risin0. Just 
as concluded in the previous works ^ the fractional 
Josephson effect occurs in the situation of even FP, 
otherwise only the normal Josephon effect can be 
observed. 

In Fig.2 we suppose Fq = Fi = 0 and choose 
Wa = 0.25 and r 2 a = 0.05 to investigate the Joseph¬ 
son effect in the case where the Majorana doublets 
couple indirectly to each other via a QD. The re¬ 
sults are shown in Figs.2-3: Fig.2 corresponds to 
the noninteracting results, and Fig.3 describes the 
influences of the intradot Goulomb interaction on 
the Josephson effects in different FPs in the case of 
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FIG. 2: Josephson current spectra in the case wht^it^ 
the Majorana doublets couple indirectly via a QD. The 
structural parameters are taken to be Wa = 0.25 and 
^ 2 a = 0.05. The left and right columns correspond to 
the even-FP and odd-FP results, respectively, (a)-(b) 
The case of the non-magnetic QD. (c)-(d) The case of 
finite magnetic field on the QD with R = 0.1. (e)-(f) 
R = 0.3. (g)-(h) R = 0.5. 



FIG. 3: Josephson current in the case where the Majo¬ 
rana doublets couple indirectly via a finite-Goulomb QD. 
The Goulomb strength is U = 2.0 and the others are the 
same as those in Fig.2. The left and right columns corre¬ 
spond to the even-FP and odd-FP results, respectively. 
(a)-(b) The case of the non-magnetic QD. (c)-(d) The 
case of finite magnetic field on the QD with R = 0.1. 
(e)-(f) R = 0.3. (g)-(h) R = 0.5. 


U = 2.0. First, in Fig.2(a)-(b) we find that when a 
non-magnetic QD is presented, it induces the occur¬ 
rence of the normal Josephson effects and the depar¬ 
ture of eo from zero weakens the current amplitudes, 
irrelevant to the FP difference. Also, the FP plays 
an important role in affecting the Josephson effects. 
Concretely, the Josephson currents in different FPs 
flow in the opposite directions for the same 0 , and 
the amplitude of is about one half of that of 

and when |5o| > 0.5 gets close to zero. Be¬ 
sides, at the points of 0 = (2n — l) 7 r, in the even-FP 
case the discontinuous change of the Josephson cur¬ 
rent is more well-defined compared with the odd-FP 
case. Next, when finite magnetic field is applied on 
the QD, the even-FP Josephson current shows little 
change except that its amplitude becomes less de¬ 
pendent on the QD-level shift. However, in the odd- 
FP case, the Josephson current changes completely. 
It can be clearly found that with the strengthening 
of the magnetic field, the original current oscillation 
is suppressed. Especially in the vicinity of 0 = 4 n 7 r, 
the current amplitude tends to disappear at the case 
of R = 0.5. Thus, it is certain that in the case of 
odd FP, a nonzero field on the QD can induce the 
occurrence of the fractional Josephson current. In 
addition, the increase of R enhances the current os¬ 
cillation around the points of 0 = ( 2 n — l) 7 r when 
5o departs from ^ 0. Up to now, we can conclude 
that when the Majorana doublets are coupled by a 
magnetic QD, the fractional Josephson effect has an 
opportunity to take place, but it is different from the 
case where the Majorana doublets couple directly.— 

Coulomb interaction is a key factor to influence 


the characteristics of QD. In Fig.3 we consider the 
case of finite intradot Coulomb interaction and in¬ 
vestigate the effect of the magnetic QD on the 
Josephson currents in the finite-Coulomb case. In 
Fig.3(a)-(b) we first find that in the even-FP case 
with a non-magnetic QD, the Coulomb interaction 
benefits the Josephson effect, since in the region of 
—2.5 < Sq < 0.5 the current amplitude is relatively 
robust and weakly dependent on the shift of the QD 
level. In contrast, for the odd-FP case, the intradot 
Coulomb interaction only moves the current max¬ 
imum to the point of Sq = — 1 . 0 , but it does not 
vary the current oscillation manner compared with 
the noninteracting case. Hence, the Coulomb in¬ 
teraction only adjusts the effect of the QD level on 
the Josephson effects but does not modify the cur¬ 
rent oscillation manner with the change of (j). Next, 
Fig.3(c)-(h) show that regardless of the FP differ¬ 
ence, the current amplitudes are suppressed by the 
application of finite magnetic field on the QD. In 
the even-FP case, the current amplitude around the 
point of eo = — 1.0 undergoes a relatively-apparent 
suppression. For the odd-FP case, except the sup¬ 
pression of the current amplitude, the fractional 
Josephson effect becomes weak but can still be ob¬ 
served. 

According to the results in the above paragraphes, 
when the Majorana doublets couple indirectly via 
a magnetic QD, the fractional Josephson effect has 
an opportunity to come into being in the odd-FP 
case. Otherwise, only the normal Josephson effect 
can be observed. In order to explain the results in 
Figs.2-3, we would like to compare the case of T 7 ^ 0 
and Va = 0 with the case of T = 0 and lU 7 ^ 0 . 
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In the former case, the Josephson effects are only 
determined by the FP of state \nfnj)f. Namely, 
when the system is located at the states |00)/ or 
|11)/, the fractional Josephson effect occurs. How¬ 
ever, when the Majorana doublets couple indirectly 
via a QD, the Fock space defined by \nfnj)f just 
becomes a subspace of the Fock space formed by 
InsTisUfTi^). Since both Fq and Fi are equal to zero 
in such a case, the nonzero contributions of |00)/ 
and 111)/ to the Josephson effect are dependent on 
their finite indirect couplings, i.e., their simultane¬ 
ous couplings to the other states. Thus, the man¬ 
ners of the indirect couplings inevitably regulate the 
properties of the Josephson effects. For instance, in 
Eq. ([T 3|) we find that in the even-FP case, the states 
that include |00)/ and |11)/ couple indirectly in the 
way of left-right symmetry. This leads to the equal 
contributions of |00)/ and |11)/ to the Josepshon 
effect, and then the normal Josephon effect takes 
place, independent of the presence of magnetic field. 
On the other hand, in the odd-FP case, Eq. m 
shows that the coupling between the states that in¬ 
clude 100)/ and |11)/ is left-right asymmetric unless 
5s 7 ^ 5s. Consequently, if finite R is considered, 
|00)/ and |11)/ make different contributions to the 
Josephson effect, which gives rise to the fractional 
Josephson effect. Also, note that in odd-EP case, 
the left-right asymmetric coupling manner weakens 
the quantum coherence and suppresses the current 
amplitude to some degree. When the Coulomb inter¬ 
action is taken into account, the QD level splits. In 
the even-EP case, the picture in the noninteracting 
case can be doubled since the left-right symmetric 
coupling manner. Nevertheless, in the odd-EP case, 
only in the case of electron-hole symmetry, the left- 
right symmetric coupling manner can be satisfied. 
Therefore, the maximal Josephson current occurs at 
the point of 5o = —1.0 in the case of U = 2.0. 

We next proceed to pay attention to the Joseph¬ 
son effect in the case where the direct and indirect 
couplings between the Majorana doublets coexist. 
The results are shown in Eigs.4-5 where Fq is taken 
to be 0.5 and Fi = 0.1. The noninteracting results 
are presented in Fig.4 and Fig.5 describes the case 
of U = 2.0. In Fig.4, we find that for any 5o, the 
two kinds of Josephson currents show dissimilar os¬ 
cillations with the adjustment of superconducting 
phase difference. In the even-FP case with = 0, 
when 5o gets approximately close to 0.25, the am¬ 
plitude of the Josephson current is decreased, other¬ 
wise, the Josephson effect will be enhanced and then 
holds. However, the current period always keeps to 
be 27r, as shown in Fig.4(a). Next, in Fig.4(b) where 
R = 0.3, we can see that the only effect of the mag¬ 
netic field in the QD is to further suppress the min¬ 
imum of the Josephson current. Such a result is 
exactly similar to the result of Fq = 0. On the other 
hand, for the odd-FP case, in Fig.4(c), it shows that 
in the region of 5i < —1.0, the Josephson current 
seems to oscillate in tt period. When the QD level 
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FIG. 4: Josephson current spectra in the case where 
the direct and indirect couplings between the Majorana 
doublets co-exist. The relevant parameters are Fq = 0.5, 
Pi = 0.1, Wa = 0.25, and r 2 a = 0.05. (a)-(b) The 

even-FP case with non-magnetic and magnetic QDs with 
R = 0.3. (c)-(d) The odd-FP results with = 0 and 
R = 0.3. 

increases to Sq = —1.0, the current period experi¬ 
ences the 7r-to-27r transition with the disappearance 
of the current around the points of (/) = 2niT. Next, 
in the region of—1.0 < 5o < 1.0, the Josephson cur¬ 
rent varies in 27r period with its maximum in the 
vicinity of 5o = 0. When sq further increases from 
1.0, the Josephson current recovers the 7r-period os¬ 
cillation gradually. Fig.4(d) presents the effect of 
the magnetic field on the QD in odd-FP case. It 
seems that in such a case, the magnetic field can not 
induce the fractional Josephson effect, but it tends 
to enhance the current amplitude in the region of 
— 1.0 < 5o < 1.0, which is exactly opposite to the 
case of Fq = 0. 

Following the above result, we present the influ¬ 
ence of the magnetic field on the finite-Coulomb 
results, as displayed in Fig.5. Here, the Coulomb 
strength is also taken to he U = 2.0. First, Fig.5(a) 
shows the even-FP result with the non-magnetic 
QD. We can find that in such a case, the current 
minimum is shifted to the position of Sq ~ —0.25. 
Besides, the Coulomb interaction efficiently weakens 
the Josephson effect, since increasing 5o from —2.0 
begins to eliminate the current amplitude gradually. 
Next when finite magnetic field is applied on the QD 
with = 0.3, it further suppresses the minimum of 
the Josephson current, similar to the noninteracting 
case [See Fig.5(b)]. The odd-FP results are shown 
in Fig.5(c)-(d), where the magnetic field strength is 
taken to be zero and 0.3, respectively. In Fig.5(c), 
we see that different from the noninteracting result, 
the 27r-period oscillation of the current occurs from 
5o = —3.0. In the region of —3.0 < 5o < 1.0, the 
Josephson current varies in 27r period with its maxi¬ 
mum in the vicinity of 5o = —1.0. Besides, it can be 
noted that the Coulomb interaction enhance the am¬ 
plitude of the Josephson current, in comparison with 
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FIG. 5: Josephson current spectra in the case of si¬ 
multaneous direct and indirect couplings between the 
Major ana doublets. The Coulomb strength is hxed at 
U = 2.0, and the other parameters are identical with 
those in Fig.4. (a)-(b) The even-FP case with non¬ 

magnetic and magnetic QDs with R = 0.3. (c)-(d) The 
odd-FP results with R = 0 and R = 0.3. 


the noninteracting current results. For the effect of 
the magnetic field in the odd-FP case, as shown in 
Fig.5(d), it is analogous to that in the noninteract¬ 
ing case. Namely, it tends to enhance the ampli¬ 
tude of the 27r-period current, despite in the region 
of—3.0 < So < 1.0. Also, it does not induce the 
appearance of the fractional Josephson effect. 

The results in Figs.4-5 can be explained similar to 
the discussion about Figs.2-3. In the case of finite T, 
the underlying physics that governs the Josephson 
effects becomes complicated. The reason arises from 
two aspects. First, the fermion number of the QD 
re-regulates the FP of \nfn^) f for conserving the FP 
in the whole system. Second, the Fano interference 
can be induced by the direct and indirect couplings 
between the Majorana doublets. Based on this idea, 
we see that in the even-FP case with zero magnetic 
field, when the QD level is far away from the en¬ 
ergy zero point, both and ris are equal to 1 or 0 
simultaneously. Consequently, the states |00)/ and 
111)/ contribute equally to the Josephson effect, so 
the normal Josephson effect occurs with the current 
amplitude proportional to Fq. Alternatively, in the 
odd-FP case with ^ 0, the 7r-period current 

occurs and its amplitude is related to Fi, due to the 
co-contribution of states |10)/ and |01)/. When the 
QD level gets close to the energy zero point, it will 
be partly-occupied. In such a situation, |10)/ and 
|01)/ contribute to the even-FP Josephson current, 
whereas |00)/ and |11)/ devote themselves to the 
odd-FP current. However, due to Fi ^ Fq, the sup¬ 
pression of only appears in a narrow region near 
the point of 5o = 0, while the 27r-period oscillation 
of distributes in a wide region accompanied by 
its enhanced amplitude [See Fig.4(a) and Fig.4(c)]. 
Next, in the presence of the intradot Coulomb inter- 




n 

r 

0 

1 2 ! 

(j)! TT 

5 4 





FIG. 6: Josephson currents in the limit of strong mag¬ 
netic held on the QD. The left and right columns de¬ 
scribe the even and odd FP results, (a)-(b) The cases of 
Fo = 0. (c)-(d) Results of Fo = 0.5. 


action, splits into two, i.e., and + U. Ac¬ 
cordingly, in the energy region of —U <63 <0, the 
fermion in the QD is changeable between 0 and 1, 
which magnihes the transformation of the Joseph¬ 
son effect caused by the shift of QD level. Since a 
hnite magnetic held on the QD plays a similar role 
in affecting the fermion occupation in the QD, the 
magnetic held makes a similar contribution to the 
Josephson effect, compared with the Coulomb inter¬ 
action. In addition, it should be noticed that the 
Fano interference induces the asymmetric spectra of 
the Josepshon currents vs 

At last, we would like to pay attention to the ex¬ 
treme case of strong magnetic held where only one 
level (i.e.. Eg) contributes to the Josephson effects. 
In such a case, the matrix dimension of and 

will be halved, as discussed in the above sec¬ 
tion. The corresponding numerical results are shown 
in Fig.6. First, in Fig.6(a)-(b) we can see that in the 
case of Fo = 0, the Josephson currents in different 
FPs are the same as each other, with their 27r pe¬ 
riod. On the other hand, when the direct coupling 
between the Majorana doublets is considered, the 
Josephson currents become dependent on the FP. As 
shown in Fig.6(c), in the even-FP case, increasing 
can change the current period from tt to 27r with the 
clear transition region near ~ —1.0. However, in 
the odd-FP case, similar result occurs when de¬ 
creases. These results can certainly be clarified by 
discussing the influence of the fermion number in the 
QD on the FP of states \nfnj)f. Also, in Fig.6(c)- 
(d) we can find that the Fano interference induces 
the dissimilar transition behaviors of the Josephson 
currents in different FPs. 


IV. SUMMARY 

In summary, we have investigated the Joseph¬ 
son effect in the Majorana-doublet-contributed junc¬ 
tion modified by the different inter-doublet coupling 
manners. It has been found that an embedded QD 
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in this junction plays a nontrivial role in modify¬ 
ing the Josephson effects, since the tunable fermion 
occupation in the QD re-regulates the FP of the 
Majorana doublets for conserving the FP in whole 
system. As a result, the 47r-period, 27r-period, and 
TT-period Josephson currents have opportunities to 
come into being, respectively. To be concrete, when 
the Majorana doublets couple indirectly via a non¬ 
magnetic QD, the normal Josephson effects occur, 
and the FP change just causes the reversal of the 
current direction and the variation of the current 
amplitude. Only in the odd-FP case, can applying 
finite magnetic field on the QD induce the appear¬ 
ance of the fractional Josephson effect. When the 
direct and indirect couplings between the Majorana 


doublets coexist, no fractional Josephson effect takes 
place, regardless of finite magnetic field on the QD. 
Moreover, the 7r-period current has an opportunity 
to appear with the shift of the QD level. We be¬ 
lieve that this work will be helpful for describing 
the QD-assisted Josephson effects between the Ma¬ 
jorana doublets. 
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